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ABSTRACT 

Planet-planet scattering is the leading mechanism to explain the large eccentricities of the observed 
exoplanet population. However, scattering has not been considered important to the production 
of pairs of planets in mean motion resonances (MMRs). We present results from a large number 
of numerical simulations of dynamical instabilities in 3-planet systems. We show that MMRs arise 
naturally in about five percent of cases. The most common resonances we populate are the 2:1 
and 3:1 MMRs, although a wide variety of MMRs can occur, including high-order MMRs (up to 
eleventh order). MMRs are generated preferentially in systems with uneven mass distributions: the 
smallest planet is typically ejected after a series of close encounters, leaving the remaining, more 
massive planets in resonance. The distribution of resonant planets is consistent with the phase-space 
density of resonant orbits, meaning that planets are randomly thrown into MMRs rather than being 
slowly pulled into them. It may be possible to distinguish between MMRs created by scattering 
vs. convergent migration in a gaseous disk by considering planetary mass ratios: resonant pairs of 
planets beyond 1 AU with more massive outer planets are likely to have formed by scattering. In 
addition, scattering may be responsible for pairs of planets in high-order MMRs (3:1 and higher) that 
are not easily populated by migration. The frequency of MMRs from scattering is comparable to the 
expected survival rate of MMRs in turbulent disks. Thus, planet-planet scattering is likely to be a 
major contributor to the population of resonant planets. 

Subject headings: planetary systems: formation — methods: n-body simulations 



1. INTRODUCTION 

The current sample of exoplanets exhibits several in- 
teresting dynamical features (Butler et al. 2006) : here 
we focus on two of these. First, there is a vast range 
of planetary eccentricities, from zero to > 0.9, with a 
median of 0.2 (0.27 for planets past 0.1 AU that have 
not been affected by tides; Rasio et al. 1996; Jackson et 
a/. 2008). Second, mean motion resonances (MMRs) in 
multiple planet systems appear to be relatively common 
(e.g., Marcy et aZ. 2001). There are 31 currently-known 
multiple planet systems comprising 44 pairs of adjacent 
planets, of which ten (23%) show some evidence of reso- 
nances (Table 1). However, the evidence for resonances 
is tentative for all but a few cases. 

Dynamical instabilities in systems of two or more plan- 
ets can explain the wide eccentricity distribution of exo- 
planets (Rasio & Ford 1996; Weidenschilling & Marzari 
1996; Lin & Ida 1997). Instabilities arise on timescales 
that are related to the planets' initial separation (Marzari 
& Weidenschilling 2002; Chatterjee et al. 2008), and lead 
to close encounters between planets and subsequent ejec- 
tions or mergers. In the aftermath of close encounters, 
the surviving planets can statistically reproduce the ob- 
served eccentricity distribution of exoplanets (Adams & 
Laughlin 2003; Juric & Tremaine 2008; Chatterjee et 
al. 2008). 
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MMRs are thought to arise primarily from conver- 
gent migration in gaseous protoplanetary disks (Snell- 
grove et a/. 2001; Lee & Peale 2002). Indeed, models 
show that capture in the 2:1 and 3:2 MMRs is a par- 
ticularly common occurrence (Thommes 2005; Pierens 
& Nelson 2008; Lee et aL2008). However, MRI-derived 
turbulence can act to remove planets from resonance. 
Adams et al. (2008) estimate that only 1% of resonant 
systems should remain for a disk lifetime of 1 Myr. 

In this paper we attempt to reconcile the planet-planet 
scattering scenario with the population of resonant ex- 
oplanets. We numerically investigate dynamical insta- 
bilities in systems of three planets located at ~ 2-10 
AU with a variety of mass distributions. We find that 
MMRs are a common occurrence, arising in 5-10 per- 
cent of unstable systems. Our simulations populate a 
range of MMRs, including the 2:1, 3:2, 3:1, 4:1 and ex- 
tending up to much higher-order (Table 2). MMRs are 
populated by scattering at random into stable regions; 
the density of resonant orbits (i.e., the fraction of phase 
space that undergoes resonant oscillations) is consistent 
with the scattered resonant systems. We propose several 
ways to discriminate between scattering and convergent 
migration as the source of exoplanet MMRs. 

2. METHODS 

Our simulations started with three planets randomly 
separated by 4-5 mutual Hill radii {RH,m = 0.5(ai -|- 
02) ([Ml -I- M2]/SMi,y/'^, where a is the semimajor axis 
and M the mass). This spacing was chosen to pro- 
duce instabilities on timescales of at least the ^ 10^ 
year timescale of runaway gas accretion^ (Pollack et 

^ Indeed, instabilities occurred on timescales from 100 years to 
98 Myr with a median of a few xlO^ years. In addition, about 1/4 
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Fig. 1. — Cumulative eccentricity distribution of the known exo- 
planets beyond 0.1 AU (thick grey), compared with our scattering 
simulations. 

al. 1996; Marzari & Weidenschilling 2002). The outer- 
most planet was placed two Hill radii interior to 10 AU; 
cases with more massive planets and therefore larger Hill 
radii therefore had the innermost planet closer to the star 
than for cases with lower-mass planets (see below) . Plan- 
ets were given zero eccentricity and mutual inclinations 
of less than 1 degree, and the stellar mass was 1 Mq. 

We considered a range of planetary mass distributions. 
For our two largest sets (1000 simulations each) we ran- 
domly selected planet masses according to the observed 
distribution of exoplanet masses: dN/dM oc M~^-^ 
(Butler et aL2006). In the "mixedl" set we restricted 
the planet mass Mp to be between a Saturn mass Msat 
and three Jupiter masses Mjup- For our "mixed2" set, 
the minimum planet mass was decreased to 10 M®. We 
also performed four "Mequal" sets (500 simulations each) 
with equal mass planets for Mp — 30 M^, Ms at, Mjup, 
and 3Mj„p. Finally, the "Mgrad" sets (250 simulations 
each) contained radial gradients in Mp. For the JSN 
set, in order of increasing orbital distance, Mp — Mj^p, 
Msat, and 30 M^. For the NSJ set, these masses were 
reversed, i.e., the Mjup planet was the most distant. The 
3JJS and SJ3J sets had, in increasing radial distance, Mp 
= 3Mj„p, Mjup and Msau and Mp = Ms at, Mjup and 
3Mj„p, respectively. 

Each simulation was integrated with the hybrid version 
of the Mercury integrator (Chambers 1999). All planets 
were assigned physical densities of 1.3 gcm~^ and col- 
lisions were treated as inelastic mergers. We used a 20 
day timestep which tests show introduces an error of less 
than 1 part in 10^ for perihelion distances larger than 0.5 
AU. In almost all cases energy was conserved to better 
than one part in 10^ for the entire 100 Myr simulation, 
which Barnes & Quinn (2004) showed is adequate pre- 
cision to test stability. However, in some cases, energy 
was poorly conserved; those cases were rerun with a 5 day 
timestep. After this step, simulations with poor energy 
conservation were removed from the analysis. 

3. RESULTS 

Figure [T] shows that four of our sets of simulations 
match the exoplanet eccentricity distribution - mixedl, 
MequaliJup, MequaliSat, and Mequal:30 M^ - with P 
values from K-S tests greater than 0.01. However, given 

of simulations were stable for 100 Myr which shows that we started 
close to the stability boundary. 
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Fig. 2. — Evolution of a system that produced a pair of planets in 
the 3:1 MMR. Top: The three planets' semimajor axes a, perihelia 
q and aphelia Q. The inner (green), middle (black), and outer (red) 
planets are 43, 105, and 16 Mgj, respectively. Bottom: Evolution 
of the 3:1 resonant argument 83 = 3A2 — Ai — (roi -|-ro2). Resonant 
libration starts immediately after ejection of the outer planet. 

the increasing number of low-mass exoplanets, we believe 
that our mixedl and mixed2 simulations are the most 
realistic initial conditions. If scattering is the source of 
exoplanet eccentricities, then soon-to-be-discovered sys- 
tems with lower-mass planets should indeed tend to have 
lower eccentricities (Ford & Rasio 2008). 

We found MMRs by examining resonant arguments for 
simulations which produced pairs of planets with period 
ratios close to commensurate values. A pair of planets is 
in resonance if any resonant argument 6i librates rather 
than circulates. For MMR p-l-q : p, arguments are of the 
form 

= {p + q)>^i - p>^2 - q^ia (1) 

where A are mean longitudes, zu are longitudes of peri- 
center, and subscripts 1 and 2 refer to the inner and outer 
planet, respectively (e.g., Murray & Dermott 1999). 

Figure [2] shows the evolution of a typical simulation 
that created a resonant system. The instability started 
55.8 Myr into the simulation, causing a series of close en- 
counters. Within a few hundred thousand years the outer 
planet was ejected and the inner two planets swapped 
places. The two remaining planets are on stable orbits 
in the 3:1 MMR, and all three resonant arguments 6*1^2,3 
librate with amplitudes between 120° and 160°. 

A variety of MMRs is populated by scattering (see Ta- 
ble 2). Most common are the 2:1 and 3:1, but higher- 
order MMRs exist up to eleventh order (13:2). The reso- 
nant libration amplitudes tend to be large, with a median 
of 110° and several cases with amplitudes of ~ 170°. This 
contrasts with MMRs from migration which are created 
in a dissipative environment and should be much smaller. 
MMRs occur preferentially in cases with mixed mass dis- 
tributions, especially those with a positive mass gradient 
such as Mgrad:NSJ. MMRs are relatively rare for equal 
mass planets, and they tend to arise more often after 
collisions rather than ejections, which contrasts with the 
mixed and Mgrad cases. This may explain why MMRs 
have not been found in previous studies (except for iso- 



Mean Motion Resonances from Planet-Planet Scattering 



TABLE 1 

Candidate Resonant Planetary systems-^ 



M(inner)/M(outer) = 
1/3 1 3/1 
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System 


0.1, a2 


ei, 62 


Ml , M2 


MMK 





(pair) 


/ A T T\ 

(AU) 




(Mjup) 






GJ 876 c-b 


0.13, 0.2078 


0.27, 0.025 


0.56, 1.935 


2:1 


II 


HD 73526 b-c 


0.66, 1.05 


0.19, 0.14 


2.9, 2.5 


2:1 




HD 82943 c-b 


0.746, 1.19 


0.359, 0.219 


2.01, 1.75 


2:1 


HD 128311 b-c 


1.099, 1.76 


0.25, 0.17 


2.18, 3.21 


2:1 




fi Arae d-b 


0.921, 1.497 


0.067, 0.128 


0.522, 1.676 


2:1 


CD 


GJ 317 b-c 


0.95, 2.35 


0.19, 0.42 


1.2, 0.83 


4:l2 


HD 108874 b-c 


1.051, 2.68 


0.07, 0.25 


1.36, 1.018 


4:1 




HD 17156 b-c 


0.159, 0.481 


0.6717, 0.136 


3.111, 0.063 


5:1 




HD 202206 b-c 


0.83, 2.55 


0.435, 0.267 


17.4, 2.44 


5:1 


in 


HD 208487 b-c 


0.49, 1.8 


0.32, 0.19 


0.45, 0.46 


7:1 



d 



"See http : / / www. Ipl . arizona.edu /~rory/research/xsp /dynamics / 1 
'Johnson ei al. (2007) did not determine apsidal angles, but 

Barnes & Greenberg (2008) used a stability analysis to predict 

that the system must be in the 4:1 MMR. 

lated cases in Adams & Laughlin 2003 and Chatterjee et 
al. 2008). 

MMRs appear to be populated at random: any stable 
region of parameter space can be accessed by scatter- 
ing. To test this hypothesis, we calculated the phase 
space density of resonant orbits within 10% of the 2:1 
and 3:1 MMRs for planetary mass ratio of 1/3, 1, and 
3, with Minner + Mouter = 400 M®. For cach MMR we 
ran ^ 22,000 3-body (star -I- two planets) simulations 
for 1 Myr. The semimajor axis of the inner planet was 
fixed at 5 AU (2:1 MMR) or 4 AU (3:1 MMR). We sam- 
pled four parameters: the orbital period ratio, the inner 
and outer planets' eccentricities, and the relative apsidal 
orientation. Inclinations (of < 1°) and mean longitudes 
were sampled at random. Resonant orbits were found by 
libration of resonant arguments. 

The density of 2:1 resonant orbits is higher for a more 
massive outer planet (Figure [3]).^ Almost all of the 2:1 
and 3:1 resonant orbits from scattering are found in ar- 
eas of high resonant density, and near-resonant "false 
alarms" lie in areas of low density. Thus, scattering does 
indeed appear to populate MMRs at random. This ex- 
plains why the 2:1 MMRs from our mixedl and mixed2 
simulations, with no initial mass gradients, have a me- 
dian Minner /Mouter of 0.5. In addition, the Mgrad:NSJ 
{Minner/ Mouter ~ 1/3) systems formed a large number 
of 2:1 MMRs while the Mgrad:JSN {Manner / Mouter ~ 3) 
cases formed far fewer. 

For the parameter space we sampled, the integrated 3:1 
resonant density is ~ 40% less than the 2:1 density. How- 
ever, the available parameter space is not evenly popu- 
lated by scattering. For example, scattering causes more 
massive planets to be closer to the star (Chatterjee et 
al. 2008). Indeed, the integrated 2:1 resonant density for 
Minner / Mouter = 1/3 is 66% higher than the 3:1 density. 
This explains the increased number of 2:1 vs. 3:1 MMRs 
in our sample - 74 cases of the 2:1 MMR and 47 of the 
3:1 (61% more in 2:1). 

The density of resonant orbits can explain other fea- 
tures of the population of resonant planets. The reso- 
nant planets tend to have smaller eccentricities than the 

" For a more detailed study of the 2:1 MMR, see Marzari et 
al. (2006) and Michtchenko et al. (2008a, 2008b). 
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Fig. 3. — Resonant density for the 2:1 MMR. For these simu- 
lations, the inner planet's semimajor axis was fixed at 5 AU, and 
its eccentricity was varied between 0.05 and 0.25. The color cor- 
responds to the fraction of orbits in resonance for a given value of 
the outer planet's semimajor axis and eccentricity (see color bar), 
averaged over 8 simulations with different apsidal alignments. 



non-resonant planets (Fig. [T|). Indeed, resonant systems 
underwent an average of about five times fewer close en- 
counters between planets before the system stabilized (~ 
30-50 vs. 100->200 encounters), compared with the me- 
dian outcome. The time between the first instability and 
stabilization for resonant cases was ^ 50,000 years, a 
factor of about five shorter than for the non-resonant 
systems. 

Low-order MMRs preferentially arise in systems con- 
taining a planet with a large Safronov number S. 
S is the ratio of the escape speed from a planet's 
surface to the escape speed from the system, S = 
{Mp/M^^y^^ {ap/RpY^^, where M^ is the stellar mass, 
Op and Rp are the planet's orbital distance and radius, 
respectively (Safronov 1969). Planets with larger S give 
stronger velocity kicks and thereby reduce the number of 
encounters needed to eject a planet. Indeed, the 2:1 and 
3:1 MMRs correlate with systems with at least one plan- 
etary S value above 4. In the mixed2 set, which is the 
only set with a significant range in 5, the median Smax 
for the 2:1 and 3:1 MMRs is 4.5, as compared with 3.5 for 
all unstable mixed2 simulations. A K-S test shows that 
the two samples are indeed different at the 99.9% confi- 
dence level. This constrains where the 2:1 or 3:1 MMRs 
can arise as a function of Mp, M^,, ap, and Rpi only sys- 
tems with relatively high-mass planets {Mp > Mjup) can 
generate these MMRs close-in. In contrast, high-order 
(4:1 and higher) MMRs from the mixed2 set match the 
sample of unstable cases and so are not constrained. 

The MMRs we found are numerically robust. The me- 
dian fractional integration error dE/E for the 170 res- 
onant systems was 2.6 x 10~^, far below the ~ 10~^ 
limit for determining stability (Barnes & Quinn 2004), 
and smaller than the median dE/E of 1.2 x 10"'' for all 
unstable systems. MMRs tend to arise in cases with 
short encounter times and relatively low final eccentric- 
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TABLE 2 

Resonances from scattering simulations 



Set 


Nsiiris 


N(%) in 


MMRs 




unstablc(%) 


MMRs 


(%) 


Mixed 1 


Qfifi-fifiq ^59%'! 


27 C4 7%") 


2-1 ("1 6%") 3-1 




fl 6%~) 4-1 
(0 7%) 5-1 
6:1, 7:2, 9:2 


Mixed2 


982-744 (76%) 


52 (7%) 


3:2' (0.8%), 2:1 




(2 4%~) 3-1 
(1.1%), 5:3, 
4:1 5:2 5:1 
7:3, 6:1, 7:2, 
8:3^ 9:4,' 10:3^ 
12:5, 11:2, 14:5 


IVLkH^IXkxL.OIVI J 


OUO ^_U^J /O ) 




7:1 


Mequal: Mj 


452-232 (51%) 


4 (1.7%) 


2:1, 3:1, 4:1 




(0.9%) 


iVl Ud.1 .IVlgfii 






Z.J. ^^±.y/o_/, O.J- 




(0.6%), 5:2, 5:1 
(C\ fi%^ 7-1 


Mequal:30Mffi 


367-365 (99%) 


10 (2.7%) 


2:1 (1.9%), 3:1, 
11:6 


Mgrad:JSN 


250-206 (82%) 


13 (6.3%) 


2:1 (1%), 3:1 




(1.9%), 4:1, 
5:2, 5:1, 7:3, 
8:3, 13:2 


Mgrad:NSJ 


245-221 (90%) 


30 (14.6%) 


2:1 (9.2%), 3:1 
(2.3%), 5:2, 
7:3, 7:2, 11:5 


Mgrad:3JJS 


250-150 (60%) 


4 (2.7%) 


3:1, 4:1 (1.3%), 
11:3 


Mgrad:SJ3J 


245-219 (89%) 


16 (6.5%) 


3:1 (5.7%), 4:1 



ities. Those situations yield smaller dE/E than for the 
more common stronger encounters that lead to very ec- 
centric planets. Thus, our cutoff of dE/E < 10~^ allows 
us to accurately sample both eccentricities and MMRs. 



4. DISCUSSION AND CONCLUSIONS 

Planet-planet scattering creates MMRs. The typical 
path to a resonant system involves several close encoun- 
ters between one smaller and two larger planets. After a 
relatively short time of instability, the smaller planet is 
destroyed, usually via ejection (78% of all cases) or col- 
lision (22%), leaving behind a pair of resonant planets. 
Relatively weak instabilities are probably very common 
in planetary systems; one may even have occurred in 
our own Solar System (Thommes et al. 1999). Nonethe- 
less, only a fraction of unstable systems produce reso- 
nant planets, typically 5-10%. These systems have large 
libration amplitudes and occupy a range of low-and high- 
order MMRs (Table 2). Most of these resonances are 
indefinitely stable; we integrated the 170 resonant sys- 
tems for an additional 1 Gyr and only 9 (5%) left the 
resonance, 4 cases leading to an additional system insta- 
bility. 

It may be possible to tell apart resonant exoplanets 

created via scattering from those created via convergent 
migration. In fact, only one resonant system appears to 
be inconsistent with a scattering origin due to its very 
low-amplitude libration (GJ 876; Marcy et aL2001). If 
two planets are trapped in the 2:1 or 3:2 MMR and the 



inner planet is the more massive, then migration can 
be stopped or even reversed (Massct & Snellgrove 2001; 
Crida & Morbidelli 2007). However, if the outer planet 
is the more massive then inward migration continues. 
In contrast, scattering produces planets in a variety of 
MMRs (including the 3:2 and 2:1) with a wide range in 
mass ratios and a preference for the outer planet to be 
more massive. Thus, scattering is likely to be responsible 
for systems past ~ 1 AU with 2:1 or 3:2 resonant plan- 
ets and a more massive outer planet. The HD 128311 
and /i Arae systems are good candidates for creation via 
scattering (Table 1; sec also Sandor & Kley 2006). 

Several extra-solar systems show tentative evidence for 
high-order MMRs - 4:1, 5:1 and oven 7:1 (Table 1; John- 
son et al. 2007; Correia et al. 2005; Gregory 2007; Short 
et aZ. 2008). No study to date has shown that migra- 
tion could capture planets in MMRs of higher order than 
2:1, although we encourage expanded studies of this pro- 
cess.'' Scattering produces a wide range of high-order 
resonances (Table 2). Thus, if the current candidate 
high-order MMR systems are confirmed (Table 1), then 
scattering is likely to be the responsible mechanism. 

Turbulence in gaseous disks may destroy MMRs, leav- 
ing perhaps only ^ 1% of planet pairs in resonance 
(Adams et aZ. 2008). This effect is stronger for higher- 
order MMRs. However, the timescale for MMR destruc- 
tion is sensitive to the strength of MRI turbulence which 
is very uncertain. In particular, if the MRI is fully or 
partially suppressed by the low ionization fraction in the 
inner protoplanetary disk (Gammie 1996) then the sur- 
vival prospects for resonant planets would be improved. 
Nonetheless, if only a few percent of systems remain in 
resonance, then scattering and migration may provide a 
comparable number of MMRs. 

Our simulations do not account for any dissipation. 
However, instabilities may occur while some gas remains 
in the disk (Moeckel et a/. 2008). In that case, MMRs 
could still result from scattering (Lee et al. 2008) and 
perhaps have smaller libration amplitudes. 

In conclusion, we have identified a new mechanism 
for the creation of exoplanet systems in MMRs. Un- 
fortunately the current data do not allow a conclusive 
determination of a resonance, let alone precise descrip- 
tions of the resonant argument oscillation. Nonetheless, 
our scattering model has important distinctions from the 
convergent migration model: high-order MMRs, large- 
amplitude resonant libration, and low-order MMRs with 
Minner / MouLer < 1- As the Orbital properties of exo- 
planets are better determined, it should be possible to 
distinguish between these scenarios. 

We thank Google for access to their machines. We are 
grateful to Greg Laughlin, Dimitri Veras, and an anony- 
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the NASA Postdoctoral Program administered by Oak 
Ridge Associated Universities through a contract with 
NASA. R.B. acknowledges support from NASA's PG&G 
grant NNG05GH65G and NASA Terrestrial Planet 
Finder Foundation Science grant 811073.02.07.01.15. 
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Highly-damped bodies can undergo resonant shepherding by shepherded planets do not survive in resonance, 

the 6:1 or even 8:1 MMRs (Raymond et ai. 2006; Mandell et 
al. 2007). However, as bodies grow the damping decreases, and 
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